An isospin dependent global nucleon-nucleus optical model 
at intermediate energies 



S. P. Weppner/'B R. B. Penney t, 1 G. W. Diffendale*, 1 and G. Vittorini^ 1 

'Natural Sciences, Eckerd College, St. Petersburg, FL, 33711 
(Dated: August 18, 2009) 

A global nucleon-nucleus optical potential for elastic scattering has been produced which replicates 
experimental data to high accuracy and compares well with other recently formulated potentials. 
The calculation that has been developed describes proton and neutron scattering from target nuclei 
ranging from carbon to nickel and is applicable for projectile energies from 30 to 160 MeV. With 
these ranges it is suitable for calculations associated with experiments performed by exotic beam 
accelerators. The potential has both real and imaginary isovector asymmetry terms to better de- 
scribe the dynamics of chains of isotopes and mirror nuclei. An analysis of the validity and strength 
of the asymmetry term is included with connections established to other optical potentials and 
charge-exchange reaction data. 

PACS numbers: 24.10.-i, 24.10.Ht, 25.40.Cm, 25.40.Dn, 25.40.Kv 



I. INTRODUCTION 



The fitting of global nucleon-nucleus optical model po- 
tentials (OMP) for elastic scattering has a venerable his- 
tory \^ I, |, II H, H, 0, I ! 0- Th e global optical 
potential of this work (WP OMP) was specifically de- 
signed with the next generation of radioactive beam ac- 
celerators in mind. It attempts to fit a lighter range of 
targets (12 < A < 60), including those far from stabil- 
ity, and it is over a limited projectile energy range (30 
MeV < E < 160 MeV). This research has produced one 
continuous global isospin dependent OMP which incor- 
porates both proton and neutron scattering and real and 
imaginary isovector asymmetry terms. Overall the po- 
tential of this work is recommended if one is interested in 
examining trends in the light and medium nuclei, specifi- 
cally isospin asymmetry dependencies, shell closure, and 
neutron excess in isotopes. Since the elastic potential is 
often used as a starting point in developing further in- 
elastic results, this global optical potential will give the 
researcher a consistent formulation over a wide variety of 
nuclei from which to generate multichannel calculations. 
To promote its use we provide an on-line optical potential 
calculator 12 1. 

Recently there has been the development of two other 
modern global OMPs for nucleon-nucleus elastic scatter- 
ing which cover an extensive projectile energy and tar- 
get range: the potentials of Koning and Delaroche (KD 
OMP) |i| and Madland (MD OMP) 0. Specifically the 
potential of Ref. [lj has set an impressive benchmark for 
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its extensive projectile energy range (1 keV to 200 MeV) 
and its accuracy of fit. We will make theoretical and 
calculatable comparisons of the new WP OMP to these 
noteworthy potentials throughout this article. 

In Sec. [Til we briefly discuss the theory of the global 
optical potential. A summary of our theoretical fitting 
procedure and the experimental reactions constraining 
the fit follows in Sec. Hill A generous amount of example 
calculations, given in Sec. lIVl compare the results of this 
potential with the two other contemporary optical poten- 
tials. This is followed, in Sec. El by a detailed analysis of 
the isovector asymmetric term identifying dramatic dif- 
ferences in the magnitude and character of the isovector 
asymmetry N ' — Z term of the three global OMPs as well 
as with standard t — p microscopic potentials. Reactions 
are isolated which would better constrain the isovector 
term in the future. We end with concluding remarks and 
future directions. 



II. THEORY 

The theoretical and computational aspects of creat- 
ing a global optical potential has been discussed in detail 
elsewhere, we only mention a few works which were found 
to be especially useful for this research. The seminal work 
of M. A. Melkanoff et. al. [TTj] discusses in thorough de- 
tail many of the hurdles that need to be overcome when 
attempting a computational fitting of nuclear scattering 
data. F. D. Beccetti and G. W. Greenless 0] created 
the first comprehensive and still viable global optical po- 
tential, we especially found poignant their discussion on 
weighting the various observables. Recently A. J. Koning 
and J. P. Delaroche have provide the nuclear physics com- 
munity with an ambitious optical model potential (KD 
OMP) [lj as previously mentioned. The work accompa- 
nying the KD OMP was beneficial for its clear discussion 
of theoretical issues and its extensive and comprehensive 
tables of available experimental data. 
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A. A new potential 



which has the following traditional format: 



Our complex phcncmological optical model potential 
contains the traditional volume (V), surface (S), and 
spin-orbit (SO) nuclear terms which are delineated us- 
ing the standard Woods-Saxon form factors 

fws{r,ni,Ai) = (l+cxp((r-7^)/A)r\ (1) 

where IZi is the radius parameter and Ai is the geometric 
diffusive parameter. The i is a placeholder index for the 
V (volume), S (surface) , and SO (spin orbit) designa- 
tions. The phenemological optical model potential takes 
the standard form: 

U(r,E,A,N, Z,V,MN) = 

-Vv (E, A,N,Z,V, MN) f ws (r, K V ,A V ) 
-iW v {E, A, N,Z,V, MN)f ws (r, K Vl A v ) 

+4A S V D (E, A)-ffws(r, K s , As) 
dr 

+iAA s W D (E,A,N, Z,V)-^fws(r,n s ,As) 

dr 

+ -V so {E 1 A,N, Z,V)^ r f ws {r,1l S o,Aso){\-<j) 
r dr 

+i-W so {E 1 A,N, Z,V)^fws{r,n S o,A S o){\ - <J) 
r dr 

+f C oui(r,n c ,A,N,Z,V) 7 (2) 

where the Vi and Wi are the real and imaginary potential 
amplitudes respectively and f CO ui is the coulomb term 



1+ V Ze 2 

f coul (r,K c ,A,N,Z,V) = — , r > ft c {3) 



feo U i(r,n c ,A,N, Z,T) = 
1 + V Ze 2 



2 211. 



2 2 
/ T \ 

~( 3 -^ y<K c . (4) 



For a neutron projectile this term is set to zero. The 
amplitudes, radii, and diffusive parameters have the fol- 
lowing dependent variables: 

• E - projectile nucleon laboratory energy in MeV 

• A - Atomic number of the target nucleus 

• N - Number of neutrons in the target nucleus 

• Z - Number of protons in the target nucleus 

• V — hi if projectile is a proton, -1 if a neutron 

• MN set to 1 if the target is traditionally singly magic 

- set to 2 if the target is traditionally doubly magic 

- otherwise set to 0. 

Explicitly the thirteen Woods-Saxon terms of this 
OMP are given using twenty-three distinct parameters 
where one Woods-Saxon potential term may have up to 
three of these parameters. The systematic polynomial 
formats of these terms are described below. The param- 
eters are purposely labeled to emphasize that each was 
fitted independently in the process of arriving at the final 
WP OMP, the techniques used will be detailed in Sec. lIIII 



Here are the polynomial forms of this optical potential beginning with the volume amplitudes: 

Vv = VV +V Vl A + V V2 A 2 + V V3 A 3 + V Vs E + V Ve E 2 + V Vl E 3 (5) 

+ V(N-Z) (v Vi0 +V Vil A + V Vi2 A 2 + V Vi3 A 3 + V Vi4 A 4 + V Vi5 E + V Vi6 E 2 ^ (6) 

+ MN (v Vm0 +V Vml A + V Vm2 A 2 + V Vm3 A 3 + V Vm5 E + V Vm6 £ 2 ) , (7) 

W v = W Va +W Vl A + Wy 2 A 2 + W v?1 A 3 + W Vs E + W Ve E 2 + Wv 7 E 3 (8) 

+ V(N-Z) (w Vi0 +W Vll A + W Vt2 A 2 + W Vi3 A 3 + W V%1 A 4 + W Vi5 E + W Via E 2 } (9) 

+ MN(w VmQ + W Vml A + W Vm2 A 2 + W Vm3 A 3 + W Vm5 E + W Vm6 E 2 ^j . (10) 

The surface amplitudes: 

V s - V So + V Sl A + V S2 A 2 + V S3 A 3 + V S ,E + V Se E 2 + V Sv E 3 , (11) 

W s = W So +W Sl A + W S2 A 2 + W S3 A 3 + W Ss E + W S<S E 2 + W Sl E 3 (12) 

+ V(N - Z)(Ws M +W Szl A + W Sz2 A 2 + W Sl3 A 3 + W St ,A 4 + W Sz5 E + W Ste E 2 y (13) 

The spin orbit amplitudes: 

Vso = V so „ + V SOl A + V S o 2 A 2 + V so ,A 3 + V so ,E + V SOe E 2 + V S o 7 E 3 (14) 

+ V(N - Z) (v SOltt + Vso n A + V SOt2 A 2 + V SOl3 A 3 + Vso^A 4 + V SOl5 E + V SOt6 E 2 ^ , (15) 
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Wso = W SOo + W SOl A + W S02 A 2 + W S o 3 A 3 + W SOr E + W S o 6 E 2 + W S o 7 E 3 (16) 

+ V(N - Z)(w SOi0 + W SOil A + W so%2 A 2 + W SOl3 A 3 + W SOt4 A 4 + W SOi5 E + W SOi6 E 2 ^ . (17) 

The volume radius and diffusive terms: 

TZy — Ry + Rv±A + Ry 2 A 2 + Ry 3 A 3 + Ry 5 E + Ry 6 E 2 + Ry 7 E 3 , (18) 

Ay = A Vo + A Vl A + A V . 2 A 2 + Ay 3 A 3 + Ay r E + Ay 6 E 2 + Ay 7 E 3 (19) 

+ V(N-Z)(A Vi0 +Ay 1 A + Ay 2 A 2 + Ay 3 A 3 +Ay 5 E + A Vt6 E 2 +Ay 7 E 3 y (20) 

The surface radius and diffusive terms: 

TZ S = R So + Rs 1 A + R S2 A 2 + R S . 3 A 3 + R Ss E + R Se E 2 + R S7 E 3 , (21) 

As = A S(> + A Sl A + A S2 A 2 + A S . 3 A 3 + A Ss E + A Se E 2 + A Sr E 3 , (22) 

the spin-orbit radius and diffusive terms: 

T^-so — Rso + Rso t A + Rso 2 A 2 + Rso 3 A 3 + Rso s E + Rso B E 2 + Rso 7 E 3 (23) 

+ V(N-Z) (i?so !0 + R SOll A + R so , 2 A 2 + R SOl3 A 3 + R S o i5 E + R S o l6 E 2 + i?so l7 ^ 3 ) , (24) 

Aso = A SOo + A SOl A + A S o 2 A 2 + A S o 3 A 3 + A S o 5 E + A SOe E 2 + A S o 7 E 3 , (25) 
and finally the the coulomb radius term: 

Tic = Rc +Rc\A + R C2 A 2 + R C . 3 A 3 + R C ,E + R Ce E 2 + R Cl E 3 (26) 

+ V{N - Z){^j) 3 (R Ci0 + R Cll A + Rc z2 A 2 + R Cl3 A 3 + R Cii A 4 + R Ci5 E + R Cl6 E 2 + R c „E 3 ^j . (27) 

I 



Parameters in this OMP which are at variance 
with the other OMPs are the real surface amplitude 
fEq. Hip (which is small), imaginary asymmetry (N — Z) 
fEqs. 19113117] . geometric asymmetry fEqs. l2QI24l2T)) , and 
the magic number dependent terms (Eqs. 171 1 0[) . The 
asymmetric parameters and their ramifications are dis- 
cussed in Sec. [V] The magic number terms attempts 
to better characterize the bonding that occurs in these 
closed shell nuclei. Our short term coulomb radius is also 
untraditionally energy dependent, the rational for this 
will be discussed in Sec. [Vj The twenty three parame- 
ters, of Eqs. ISTI2T1 which describe the thirteen potential 
terms of Eq. [21 are listed in Table Q] 

This potential was put into a standard optical poten- 
tial calculator which solves the Schrodinger equation for 
spin i-spin scattering using a distorted Born wave ap- 
proximation (DWBA) in a coulomb wave function ba- 
sis. TALYS [llEl was used which applies ECIS [H to 
calculate the solution once the final product was devel- 
oped. We have produced a Java applet [l2| which con- 
tains our own optical potential calculator as well as some 
useful input files for use in TALYS and ECIS which will 
let researchers produce results quickly. The gestation of 



the parameters will be described in Sec. |TTT] but first an 
overview will be given of the theoretical structure of this 
and the other recent global OMPs. 



B. Theoretical comparison with other global 
optical potentials 

The general design of all three OMPs under study 
(KD 0, MD 0], WP) is the same. They use similar 
Woods-Saxon functional forms which include volume and 
spin-orbit terms. Additionally the KD and WP opti- 
cal potentials include surface terms. The KD and MD 
potentials also include an additional coulomb correction 
term, in the WP optical potential the short term coulomb 
potential includes energy dependence to include these 
coulomb correction effects (this will be discussed in great 
detail in Sec. IV Bp . Both the WP and MD produce one 
potential that is utilized for both proton and neutron 
scattering while the KD potential has separate potentials 
for the isospin dependence of the projectile. 

The KP OMP [l[ imaginary parameters are deter- 
mined by a dispersive relationship which is dependant 
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TABLE I: Model parameters for the subject of this work, the WP global optical potential Each term is a 5 to 7 term separable 
polynomial in A are E which are given in Eqs. 151271 Tools have been developed to facilitate the use of this potential including 
an on-line optical potential calculator [l2| ]. 



on the difference between the projectile laboratory en- 
ergy and the Fermi energy of the target potential as well 
as constants designated as part of the real potential. Of 
the three optical potentials discussed it has the deep- 
est theoretical underpinnings and given these dispersive 
constraints and high accuracy it has successfully pushed 
the theoretical development of the global optical model 
potential to a new level. 

A goal of the MD OMP 0] is to describe the elas- 
tic scattering data sets using fewer parameters. It uses 
constants, linear, and the occasional quadratic forms to 
describe all of its Woods-Saxon parameters. It is quite 
impressive that a high quality isospin dependent poten- 
tial was produced with so few parameters. 

In contrast this work (WP OMP) uses quadratic, cu- 
bic, and occasionally quartic polynomials that have no 
direct relationship to formal scattering theory. It has 
the highest amount of adjustable parameters of the three 
potentials. What it achieves is ease of use, good asymme- 
try and mirror nuclei analysis, and complete separabil- 
ity between the nuclcon (A) and energy (E) parameters. 
This WP potential also includes a direct imaginary vec- 
tor isospin asymmetry (N—Z) terms which neither of the 
other potentials have. In Sec.|V]a detailed comparison of 
the isovector differences of these three optical potentials 
will ensue. 



III. PROCEDURE 

This research tried ambitiously to minimize a \ 2 (pro- 
portional to the square of the difference between the theo- 
retical fit and the experimental data) . on over three hun- 
dred different nucleon-nucleus experiments by adjusting 
the polynomial fit to the twenty three parameters given 
in Table Q] To attempt to decouple the terms from each 
other, the twenty-three values were varied using a sys- 
tematic method detailed below. A listing and discussion 
of the experimental dataset used to constrict the vari- 
ables then follows. 



A. Calculation Techniques 

The fitting elastic nucleon-nucleus scattering code was 
developed by two of the authors (S. Weppner and R. Pen- 
ney). A numerov routine found in Ref. [11] was used 
to solve the non-relativistic position space Schrodinger 
equation in a coulomb basis with a relativistic correc- 
tion found in Ref. The routine which produced the 
coulomb wave functions was found in Ref. 161 . A Powell 
routine, adopted from Numerical Recipes [171] ■ was used 
to minimize a weighted x 2 . 

The code was developed to run on a multi-processor 
parallel system. Each processor was assigned at least one 
nucleon-nucleus experimental data set at a given energy 
and nucleon number (A). This data set could be as sim- 
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pie as one experiment or many experiments including the 
observables of total neutron cross section, total reaction 
cross section, differential cross section, and polarization. 
Each set contained, if available, both proton and neutron 
observables and varying target proton numbers that all 
shared a common target nucleon number and projectile 
energy. For example A = 40, E — 40 MeV experimen- 
tal data exists for a proton striking 40 Ca and producing 
a total reaction cross section observable, a differential 
cross section, and polarization. These were fit simulta- 
neously with data that exists for a neutron striking Ca 
and producing a differential cross section, a total reac- 
tion cross section, a total cross section, and finally also 
including a proton at 40 MeV striking 40 Ar and produc- 
ing a reaction and differential cross section. All eight of 
these experiments made a complete working data set in 
which the parameters were varied and the minimization 
routine for the weighted \ 2 was executed. By analyzing 
many different data sets together at the same energy and 
nucleon number it was intended to reduce systematic er- 
ror and ultimately derive a better global fit among chains 
of target nuclei. 

Each fitting cycle was comprised of a three step pro- 
cess. First all parameters were adjusted except three 
(the two magic number and the asymmetry coulomb pa- 
rameters: Eqs. I7I10I27I were held fixed). Since the other 
twenty parameters were being adjusted simultaneously 
it was important to analyze only a subset of the entire 
experimental dataset to avoid ambiguities within the pa- 
rameter space. The datasets used encompassed either 
only N = Z targets or certain sets in which a variety 
of different Z targets existed for a fixed nucleon num- 
ber and fixed projectile energy or certain sets in which 
both neutron and proton projectile observables existed 
for that fixed A and fixed projectile energy E. These 
sets therefore either had zero N — Z dependence or this 
dependence was clearly delineated by including a variety 
of elements as targets and/or projectiles. There were 115 
data sets of fixed E and A experiments on which these 
twenty parameters were allowed to vary while searching 
for a minimum weighted x 2 . The best values for the pa- 
rameters were then fit to the polynomials of A and E 
listed in Eqs. EE3 

Second all data sets with non-zero (N — Z) depen- 
dence were used to adjust thirteen parameters. There 
were thus ten parameters held fixed which were the non- 
asymmetric potential parameters that had other asym- 
metric terms within the same Woods-Saxon potential 
term (Eqs. I5I8I12I14I16I19I23I26I and the magic number 
terms Eqs. 171101 were held fixed), the other thirteen terms 
were varied. There were 90 experiment data sets (each 
data set might contain more than one experiment as de- 
tailed above) with some N — Z dependence that were 
used for this task. The best values for the parameters 
were then fit to the polynomials in A and E. The pa- 
rameters that were free to be modified in both of these 
steps were deduced using an average of both values which 
produced the minimum weighted x 2 results. 



Lastly, the magic number terms were adjusted 
f Eqs. [TUTU were varied) to find the minimum weighted x 2 
keeping all other parameters fixed. For the magic num- 
bers the traditional 2,8,20, and 28 were used. There were 
40 data sets which contained at least one target magic 
nuclei experiment which were used to constrain these last 
two parameters. 

The variances in the parameters were minimized to a 
weighted x 2 which favored forward angles over backward 
angles (in the x calculation there was a square root re- 
lationship such that the fiftieth forward angle point in a 
differential cross section was weighted only about one sev- 
enth that of the extreme forward point following Ref. [3(). 
Differential cross sections were favored over polarizations 
by a factor of 1.5 and in general neutron total cross sec- 
tion point were favored to be approximately equal to half 
of a complete typical differential cross section. Used, but 
favored at half the weighting of the total neutron cross 
section, were the neutron and proton reaction cross sec- 
tions. Choosing the correct weightings was found to be 
an art form where a balance tenuously existed in which 
every reaction was regarded but certain reactions were 
strengthened so that the x 2 parameter space contained 
large relative minima which the search functions could 
find easily. 

To find these minima a Monte Carlo preliminary gross 
search was done using a Sobol [18| number generator. 
Since the parameter space was often twenty dimensions 
this might include up to 3 x 10 5 vectors in which weighted 
X 2 were first calculated by solving the Schrodinger equa- 
tion for the given potential. The lowest ten vectors were 
analyzed in more fine detail by seeking the local minima 
in their vicinity within the x 2 parameter space. The x 2 
minimization program sought the steepest derivative in 
the multi-dimensional parameter space [l7|. Although 
the overall quality of the fits were examined occasionally, 
the entire process was close to automatic. 

At first the twenty three parameters had large vari- 
ances which were slowly reduced following the three 
steps described as the weighted normalized per point chi- 
squared was reduced and approached a global minimum 
gracefully. The highest 5% of the x 2 were thrown out af- 
ter each fit (which routinely were the same sets) and then 
each parameter was fit originally to a cubic polynomial 
in A (nucleon number) and E (projectile energy). Even- 
tually many of the cubics were reduced to quadratics if 
warranted and in the asymmetry terms the A polynomial 
was increased to fourth order to give a better fit over the 
entire nucleon target atomic number range. 



B. Experimental Datasets 

This work strived for a comprehensive collection of ex- 
perimental data sets. As a starting point the excellent 
summary of elastic data found in Ref. [l[ was used, the 
product of the KD optical potential group. This was 
supplemented with new data and also nuclei outside the 
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Proton-Nucleus Experimental Elastic and Total Reaction Data References 


Nucl. 


[Ref] (Energies (MeV); 
Observables) 


Nucl. 


[Ref] (Energies (MeV); 
Observables) 


Nucl. 


[Ref] (Energies (MeV); 
Observables) 








20](20-84;dcs,pol) 






12 Be 


19](55;dcs) 


12 C 


|2l|(30;dcs,pol) 


13 C 


[H[(30;dcs,pol) IH](35;dcs) 


14 N 


23](30;dcs) 22](35;dcs) 
[H,ii|(50;dcs,pol) [H(35;rcs) 
H2(122;dcs) [II](142;dcs) 




[Hl](30-60;rcs) [13](35;dcs) 
[H(35;dcs) J2§](40;dcs) 
[H(49;dcs) [13(50;dcs) 




[H](35;dcs,pol) 
[23,|3l(72;dcs,pol) 

36](135;dcs,pol) 


« N 


22](35;dcs) H3](39,44;dcs) 
4l](42,44;dcs) 31] (65;dcs,pol) 




[H(61;dcs) [H(65;dcs,pol) 
g3(65;rcs) [H(75,150;dcs,pol) 


16 o 


21] (30;dcs,pol) 40] (35;dcs,pol) 
IH, H3](35;dcs) [13(30-47;rcs) 


17 


22](35;dcs) 45](66;dcs) 




[H(81-180;rcs) g3](96,dcs) 




5i](43,46;dcs) 


18 


22](35;dcs) 41] (42,44;dcs) 
:50](43;dcs) 45](67;dcs) 




[H(122,160;dcs,pol) 
[5l|] ( 1 3 5 ; dcs ,p ol , res) 




[H,|3i(49;dcs,pol) 
H2(65;dcs,pol) |i|(65;rcs) 


20 q 


53](30;dcs) 50](43;dcs) 




[H(142;dcs) Hi(145;dcs,pol) 




56](135;dcs,pol) 


20 Ne 


;22](35;dcs) ]27](35;rcs) 




[3l|(150;dcs pol) 


22q 


57](47;dcs) 




52](65;dcs,pol) 




[i|(156;dcs,pol) 


22 Ne 


22, 37](35;dcs) 


24 Mg 


[M|(30-45;dcs) [60](50;dcs,pol) 


28 Si* 


61](52;dcs) 


29 Si 


30](72;dcs,pol) 




l2i,[53(65;dcs,pol) 


30 Si 


62](52;dcs) 


Sip 


30](72;dcs,pol) 




63] (80;dcs,pol) 64] (135;dcs,pol) 


32 Si 


65](42;dcs) 


32g 


66](53;dcs) 39] (65;dcs,pol) 




67](30;dcs) 


34 Ar 


66](47;dcs) 


36 Ar 


68](33;dcs) 


37 C1 


22, 37](35;dcs) 


39 K 


22, 37](35;dcs) 


40 Ar 


69](30,33,37,41;dcs,pol) 


40 Ca* 


46](80-180;rcs) 70] (30-48;dcs) 


42 Ar 


68](33;dcs) 




H2, IH, H3 (30-50;dcs,pol) 


42 Ca 


[73(30-48;dcs) [n|](30-48;rcs) 


44 Ar 


68](33;dcs) 




27](36-47;rcs) 52](65;dcs,pol) 




H2](35;dcs) Jz|(49;dcs) 


44 Ca 


[73](30-48;dcs) [zl(30-48;rcs) 


45 Sc 


37](35;dcs) 72] (50;dcs,pol) 




[z3(65;dcs,pol) 




[72|(49;dcs) [H,[3|(65;dcs,pol) 


46 T ; 


52, 73](65;dcs,pol) ;75](100;dcs) 


48 Ca 


[73|(30-48;dcs) [zl(30-48;rcs) 


48 Ti 


76](30-48;rcs) 


50 Ti 


73](65;dcs,pol) 75] (100;dcs) 




H2,[7i,[73(65;dcs,pol) 




[H,0(65;dc S ,pol) 75 M ()():. I.'m 


50 Cr 


61](52;dcs) 73] 


52 Cr 


73](65;dcs;pol) 




[Z3(100;rcs) 


54 Cr 


22](35;dcs) 73] (65;dcs;pol) 


54 Fe * 


@(30,40,62;dcs) 




58](156;dcs) 


56 Fe* 


_24](40,60;rcs) 79](65;dcs,pol) 




[§(J(30-48;rcs) [8l|](40;dcs,pol) 


58 Fe 


22](35;dcs) 


59 Co 


82](30;dcs) [§3|(40;dcs) 
|3](30,40;rcs) 0(65;dcs,pol) 



TABLE II: This is a partial listing of the proton-nucleus experimental data used in the fitting process of this work, the remainder 
is located in Table 7 of Ref. [J. The references are listed followed by the laboratory energy of the projectile and the observables 
found in the reference article (res - total reaction cross section, dcs-differential cross section, pol-analyzing power). The data 
has an energy range from 30 MeV to 160 MeV, and a nucleon number, 12 < A < 60 as required by the calculation. Table 7 of 
Ref. [l[ contains* a substantial listing of experimental data references primarily for the 27 A1, 28 Si, 40 Ca, 54 Fe, 56 Fe, 58 Ni and 
60 Ni, all those experimental references were also used in this work if they were within this optical potentials applicable energy 
range. 



Neutron-Nucleus Experimental Elastic and Total Reaction Data References 


Nucl. 


[Ref] (Energies (MeV); 
Observables) 


Nucl. 


[Ref] (Energies (MeV); 
Observables) 


Nucl. 


[Ref] (Energies; (MeV) 
Observables) 


12 c 


85](30-49;rcs) 86](35;dcs) 


14 N 


85](30-49;rcs) 


16 


85](30-49:rcs) 56](35;dcs) 




H3](55-75;dcs) [H(65-156;dcs) 


27 Ar 


85](30-49;rcs) 


28 Si* 


1891] (misc. ;rcs) 




[§(|(96;dcs) [89l,[9H(misc.;rcs) 


40 Ca* 


]85](30-49;rcs) 88] (65-156;dcs) 


56 Fe* 


1891] (misc. :rcs) 



TABLE III: This is a partial listing of the neutron-nucleus experimental data used in the fitting process of this work, the 
remainder is located in Table 1 of Ref. The references are listed followed by the laboratory energy of the projectile and 
the observables found in the reference article (res - total reaction cross section, dcs-differential cross section). The data has 
an energy range from 30 MeV to 160 MeV, and a nucleoli number, 12 < A < 60 as required by the calculation. Table 1 of 

* 24 27 \ 28 32 40 56 

Ref. [1[ contains a substantial listing of experimental data references primarily for Mg, Al, Si, S, Ca, and Fe, all 
those experimental references were also used in this work if they were within this optical potentials applicable energy range. 
The neutron total cross section data sets of Refs. [92I . [93| were also extensively used. 
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assumed range of the KD potential. These additional 
data sets were either nuclei lighter than 27 A1 or those 
nuclei which are non-spherical and were not considered 
in that article, there were also a few additional data sets 
that were discovered and added to our database. All 
proton- nucleus data sets used that were not in Ref. [l[ 
are listed in Table HT1 The smaller set of neutron- nucleus 
data sets used here and not cited in Ref. |l[ are listed in 
TableHIIl A reader of this work, Ref. [J, and the growing 
EXFOR/CSISRS database at the National Nuclear Data 
Center [94j have a near complete compilation of elastic 
nucleon-nucleus and total cross section experimental data 
listings at intermediate energies. 

Not all the experimental data was used to help con- 
strain this OMP. Some data, usually published before 
1960, had large systematic differences with later reactions 
and were disregarded. Also, as mentioned in Sec. IIII1 
calculations were only done if angular differential data 
existed at a given energy and nucleon number. It was 
established for this fitting procedure that total cross sec- 
tion data did not have enough information (only one data 
point) to constrain the parameters adequately so it was 
only used in conjunction with differential experimental 
data. Occasionally if the energy of a total cross section 
was close to an energy where differential data existed, 
the total reaction cross section data was adjusted follow- 
ing the forms given in Refs. (95ll9^|. This procedure was 
not needed for the total neutron cross section data of 
Refs. [92|, [93[ since the energy coverage was substantial 
for these datasets. 



IV. RESULTS 

To produce the calculations shown in this section the 
DWB A scattering code TALYS [H, G3 was used for con- 
sistency. All experimental data is shown as black cir- 
cles, the new results of this work are depicted as light 
solid green lines (WP OMP), the calculations of Kon- 
ing and Delaroche 1] (KD OMP) are shown using dark 
blue dashed lines and the calculations of Madland [2| 
(MD OMP) are depicted using a medium red dot-dashed 
lines. The motivation is not to be exhaustive but to give a 
fair representative overview of the features of these three 
modern optical potentials. As a global summary it can 
be concluded that that all three potentials do a fairly 
good job fitting the presently available elastic scatter- 
ing data and dramatic contrasts are not proffered until 
Sec. El The details of this section do illustrate some mi- 
nor differences and so first their will be an examination of 
the neutron-nucleus observables then following with the 
proton-nucleus observables. 



A. Neutron-Nucleus Observables 

In Figs. 1-4 the calculations of total neutron cross sec- 
tion results to the comprehensive data sets of Refs. [92|, 



[93j are compared. Overall the calculations of the KD 
OMP of Ref. [H do the best job at reproducing the ex- 
perimental data, usually within a remarkable 1%. The 
potential of this work (WP OMP) is usually within 5% 
of the experimental data sets and that of the MD OMP 
is usually within 10%. These are plotted on an elongated 
linear scale to accentuate the disparity but it should be 
recognized that even a 10% difference between theory and 
experiment is extraordinary for most nuclear scattering 
observables. 

The remarkable fit by the KD OMP calculation is jus- 
tified, the authors considered this data important and 
thus weighted it accordingly. They also had, as discussed 
in Sec. UH a separate proton nucleus and neutron nu- 
cleus potential. At energies greater than 50 MeV, where 
neutron-nucleus elastic scattering data is scarce, this was 
one of the few observables used to fit their neutron optical 
potential. With fewer constraints this observable, with 
its high caliber of data [92|, [93[, was easier to fit. Con- 
versely the optical potential of this work and the optical 
potential of the MD OMP had to continually compro- 
mise by both fitting neutron-nucleus and proton-nucleus 
observables simultaneously. 

In Fig.[T]the calculations are fit to four standard nuclei. 
As in all these calculations the WP calculation of this 
work (solid green line) was not tested below 30 MeV and 
the MD calculation was not used below 50 MeV. Likewise 
the KD calculation is not applicable to the lighter nuclei 
(A < 27) targets. Systematic trends emerge where the 
KD calculation runs lower and closer to the experimental 
results than the other two calculations for the neutron- 
nucleus total cross section. The WP calculation is also 
systematically lower and closer than the MD calculation 
for the neutron total cross section observables as typified 
in this figure. There is a small kink in the MD results at 
E=130 MeV, this is a real discontinuity in the Woods- 
Saxon function form parameters for this potential at both 
130 MeV and 140 MeV. The history of the phenemologi- 
cal fitting endeavor shows that researchers have struggled 
to fit these higher energies well, the present works are no 
exception. 

Fig. [3 contains the lighter targets which follow the 
same trends as described for Fig. [T] Of the two appli- 
cable optical potential calculations (the KD OMP is de- 
fined only when A > 24) the WP OMP of this work does 
better in all cases except for neutron scattering from 19 F. 
The odd nuclei have non-zero spin forces which are not 
included in the solution technique which in part may ex- 
plain the anomaly. 

Concluding the results of neutron-nucleus total cross 
section calculations are Fig. 13141 which contain a variety 
of less common nuclei and odd spin nuclei neutron total 
cross sections. The KD optical potential of Ref. [l] which 
was explicitly not fit to these non-spherical nuclei still 
does remarkably well, the same trends discussed in the 
earlier figures still hold. What becomes apparent in the 
heavier nuclei (50 < A < 60) calculations of this work 
(WP OMP - solid light green line) is that there seems 
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3000 



Experiment 

KDOMP 
MDOMP 
WPOMP 



2500 



2000 ' E 



1500 



1000- 



500- 




E lab (MeV) 

FIG. 1: (Color online) The experimental total cross section for neutrons scattering from 27 Al[§i], 28 Si[92], 40 Ca[93|] and 56 Fe[93| 
from 20 MeV to 160 MeV for the laboratory energy of the neutron. They are fit to three different optical potential calculations 
(KD[l[,MD@] and to this work: WP) which are described in the legend. The experimental data and theory have been offset 
by constant amounts for multiple comparisons on one graph. The WP and MD calculations have a minimum energy limit of 
30 and 50 MeV respectively. 
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FIG. 2: (Color online) The experimental total cross section for neutrons scattering from the light nuclei of 
12 C[93, 14 N[9i], 16 0[93] and 19 F[§3[ from 20 MeV to 160 MeV for the laboratory energy of the neutron. They are fit to 
two different optical potential calculations (MD0] and to this work: WP) which are described in the legend. The experimental 
data and theory have been offset by constant amounts for multiple comparisons on one graph. The WP and MD calculations 
have a minimum energy limit of 30 and 50 MeV respectively. The KD calculation is not applicable for these light nuclei. 
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lab(MeV) 

FIG. 3: (Color online) The experimental total cross section for neutrons scattering from ^PlIl^Stlll^KjIII and 48 Ti[H|] 
from 20 MeV to 160 MeV for the laboratory energy of the neutron. They are fit to three different optical potential calculations 
(KD[l[,MD@] and to this work: WP) which are described in the legend. The experimental data and theory have been offset 
by constant amounts for multiple comparisons on one graph. The WP and MD calculations have a minimum energy limit of 
30 and 50 MeV respectively. 
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Experiment 
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MDOMP 
WPOMD 



3000- 



2500- 



2000 



1500- 



1000- 




E kb (MeV) 

FIG. 4: (Color online) The experimental total cross section for neutrons scattering from 51 V[§3, 52 Cr[9i|, 55 Mn[93] and 59 Co[93l] 
from 20 MeV to 160 MeV for the laboratory energy of the neutron. They are fit to three different optical potential calculations 
(KD[l|,MD|2j and to this work: WP) which are described in the legend. The experimental data and theory have been offset 
by constant amounts for multiple comparisons on one graph. The WP and MD calculations have a minimum energy limit of 
30 and 50 MeV respectively. 
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(deg) (deg) (deg) 

cm. e/ cm. c/ c.m. v a ' 

FIG. 5: (Color online]Neutron-nucleus elastic differential cross section experimental data for the target nuclei of 28 Si[87l.l97l.l98l|. 
40 Ca[H,|9i,llil, 56 Fe [87U97I. |9^. [lOOl] at a variety of incident laboratory energies. They are fit to three different optical potential 
calculations (KD[l|,MD@jand to this work: WP). As in all figures, KD OMP is a blue dashed line, MD OMP is a red dot-dashed 
line and WP OMP is a green solid line. The experimental data and theory have been offset by constant amounts for multiple 
comparisons on one graph. The target labels at the top read left to right correspond to the calculations and experimental data 
read top to bottom. The MD calculations has a minimum energy limit of 50 MeV and is therefore missing from calculations 
below that energy. 



to be a systematic energy shift; the shape of the curve 
is good, but it seems to contain a small shift towards 
higher energies, a possible explanation will be proffered 
in Sec. IVBl 

The neutron-nucleus differential experimental data is 
scarce above 30 MeV (see Table [Til] and Ref. [J). A rep- 
resentative sample of the data with calculations is shown 
in Fig. [SI Overall all three potentials describe the ex- 
perimental data adequately, most impressive are the re- 
sults of the calculation of the MD potential which simul- 
taneously fits both proton-nucleus and neutron-nucleus 
data while also using the fewest terms and parameters 
of the three optical potentials examined. Systematically 
the WP potential of this work has the weakest results, 
especially when the scattering angle is greater than 45 de- 
grees This implies that to improve these results the WP 
potential would need to give more weighting to higher 
angle results than was determined (see Sec. IIIII which 
discusses the weightings chosen) . To adequately measure 
the effects and form of the isospin dependent and asym- 
metry terms in the optical potential which is used by the 



WP and MD OMPs more high energy neutr on- nucleu s 
differential and reaction data is sorely needed [lOll Il02l | . 

Overall the best optical potential calculation for 
neutron-nucleus scattering is the work of Koning and Dc- 
laroche (KD OMP) [j. It has the best fit to the total 
cross section observables (even the non-spherical and odd 
nuclei which were not used to constrain that potential) 
and it also does an admirable job with the differential ob- 
servables. It has a wide energy range and is thus suitable 
for systematic neutron- nucleus studies. Its largest deficit 
is that it is not applicable to light nuclei (A < 24) which 
are important in astrophysics and biological physics. 

B. Proton-Nucleus Scattering Observables 

The proton-nucleus elastic scattering observables will 
now be considered. The potential used by the MD and 
WP calculations remains the same while the KD po- 
tential uses a different optical potential to calculate the 
proton-nucleus observables. 



13 



2500- 



2000- 



P6 



1500- 




E lab (MeV) 



FIG. 6: (Color online) Proton-nucleus total inelastic cross section data for the targets 12 C, 16 , 27 Al, 40 Ca, 56 Fe, and 58 Ni from 
20 MeV to 160 MeV for the laboratory energy of the projectile proton. The experimental data comes from a variety of sources 
and are compiled and discussed in Refs. [H. Y7n. l95j . There are also new higher energy measurements for 12 C and 40 Ca found in 
Ref. [4^]. They are fit to three different optical potential calculations (KD|i[,MD[2] and to this work: WP) which are described 
in the legend. The experimental data and theory have been offset by constant amounts for multiple comparisons on one graph. 
The WP and MD calculations have a minimum energy limit of 30 and 50 MeV respectively. 
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(deg) (deg) (deg) 

cm. 07 cm. 07 cm. 07 

FIG. 7: (Color online) P roton-nuc leus elast ic R utherford reduced differential cross section experimental data for the target 
nuclei of 27 Al[13, H \M DM Gil], 28 Si[H, Efll, 4 °Ca[|, H [nj, at varying proton laboratory energies. Refer to Fig. [5] for 
details of the legend for the theoretical calculations. 



A representative sample of the inelastic or reaction 
total cross section data are shown in Fig. [5] with pro- 
jectile energies ranging from 30 MeV to 160 MeV. The 
calculations in this case are not as close to the exper- 
imental data, and the data itself is much sparser and 
more unsure than the earlier neutron total cross section 
data. More accurate high energy data is needed in this 
observable to better constrain the optical model theory 
(it plays a s ignif icant role in determining the absorptive 
strength [ifjl. IllOj ]). an example will be proffered to illus- 
trate this in Sec. [V] All three global optical potentials 
do well but are not excellent. As with the neutron total 
cross sections the legend is the same and the figure is 
elongated and plotted on a linear scale to emphasis the 
differences. 

There is a much more substantial amount of experi- 
mental data for the proton-nucleus differential cross sec- 
tion and the spin observable analyzing power and a rep- 
resentative sample is contained in Figs. [7lTl3l In Fig. [7] 
there is a comparison of some of the common nuclei tar- 
gets, 27 A1, 28 Si, and 40 Ca for the differential cross section 
normalized to the coulomb Rutherford differential cross 
section. Then an examination of the heavier common 
target subset ( 54 ' 56 Fe and 58 N) transpires in Fig. [5J All 
these nuclei were calculated in Ref. |l| by the KD calcu- 
lation and in this work are compared with the MD and 



WP potentials (in the applicable energy ranges, an MD 
calculation is not created for an energy of the projectile 
of less than 50 MeV). Overall all the three calculations 
do well. Some systematic trends become apparent: the 
WP OMP struggles at the higher angles to reproduce the 
experimental data (as it did with the neutron-nucleus 
observables) . The disappointing fit at larger angles in- 
sinuates that the balance of the weighting functions for 
the WP OMP favored the forward angles too much. The 
KD OMP often exaggerates the minimum which prob- 
ably signifies that the coulomb strength is slightly weak 
short range as prescribed in that potential (in general the 
long term coulomb force mildly obscures the diffraction 
effect generated by the short term forces) . Because these 
are plotted on a logarithmic scale the percent difference 
between the theory and experiment is often larger than 
it appears (sometimes as much as 50% compared to 10% 
for the total cross sections). The quality of fit, as as- 
certained by an analysis of the \ 2 i n 

Ref. Q, is not as 

good as with the neutron observables, yet it is still quite 
impressive for all three optical potentials. 

The next two figures (Fig. [JJ] and Fig. [TTJ]) examine the 
reduced differential cross section of target nuclei which 
are non-spherical and farther from the line of stability 
than fit in Ref. [1J using the KD potential. Many of 
these examples are odd non-spin zero nuclei. Overall the 
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cm. 0/ cm. 0/ cm. 07 

FIG. 8: (Color online) Proton-nucleus elastic Rutherfor d reduced diff erential cross section experimental data for the target 
nuclei of 54 Fe 0,51 El El, 56Fe [l [H. H. IH. 58 NiH3, Eol, [Io3, Eol] , at varying proton laboratory energies. Refer to Fig. [5] 
for details of the legend for the theoretical calculations. 



three potentials do surprisingly well with the same sys- 
tematic problem as with the more popular and common 
nuclear targets, however some aspects are dissatisfying. 
In Fig. [9] the results of the proton- nucleus reduced differ- 
ential cross section of the calcium and chromium isotopes 
are shown and none of the calculations do extremely well 
describing all four isotopes simultaneously. The research 
of this work focused on trying to find the strength of 
the isovector antisymmetry term by simultaneously fit- 
ting these isotope's observables together. The resulting 
calculation (the solid green line) fit 40 Ca and 44 Ca iso- 
topes comfortably but struggle with the the two other 
isotopes, the other calculations are in similar predica- 
ments. The general shape is correct but the minima and 
maxima are often missed by over 25%. It seems that a 
asymmetry term which is linear and has mirror symme- 
try evades accurate discovery. The same difficulties with 
describi ng th e calcium isotopes were recently discussed 



in Ref. [10l|. These same arguments can be made, to 



a lesser extent, with chromium as shown in Fig. [5J The 
failure to create an excellent isovector asymmetry term is 
the motivation for the analysis of Sec. [V] It also must be 
recognized that these potentials are often also missing a 
general spin-spin term which may be important for non- 
spin zero targets. As this potential was being developed 
a spin-spin term was used but the results had a substan- 



tial amount of noise and it was eventually removed. A 
competent analysis of the validity and strength of this 
spin-spin term would be a worthy endeavor for future 
work. 

Fig [TU] contains the reduced differential cross section 
of some lighter targets (if A < 27 the KD calculation is 
not applicable). This figure shows the struggles the cal- 
culations have to fit at higher energies which has been a 
common problem with an optical model calculation. [49| . 
Likewise the back angles continually cause difficulty in 
the lighter targets. Recently there has been success in 
fitting these back angles by adding term s which add 
non-local approaches to these local OMPs jllll . Ill2| by 
adding angular momentum and parity dependencies to 
the potential. 

In the center panel are highlighted some newer dif- 
ferential cross section data of the oxygen isotopes. The 
WP calculation (which is the only one applicable) does 
well in describing the general trend of the first diffrac- 
tion minimum shifting as the neutron number increases. 
The mirror nuclei of calcium and argon are the targets 
in the right most panel. A recurring disappointment in 
not being able to fit all the isotopes with the same ex- 
cellent quality of fit, an important motivation for this 
work. However, for those working with exotic beams, 
this potential does give a good starting basis for contin- 
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cm. e/ cm. e/ c.m. v e/ 

FIG. 9: (Color o nline ) Pro ton- nucleus elastic Rutherford reduced differential cross section experimental data for the target 
nuclei of 34 Ar[g||, GUI, 30 S[lOl, 40 Ar[z3l, 45 Sr[z2, 40 ' 42 ' 44 ' 48 Ca[7l], 50 ' 52 ' 54 Cr[7l]], at varying proton laboratory energies. 
The target labels at the top read left to right correspond to the calculations and experimental data read top to bottom for the 
two rightmost panels. Refer to Fig. [S] for details of the legend for the theoretical calculations. 



uing research. 

In Figs. [TTIfT^l are plotted a representative sample of 
the proton-nucleus spin analyzing power (polarization) 
observables using the same legend as in all earlier figures. 
In Fig. [IT] which examines tradition nuclear targets, all 
the calculations do well. The fitting of the polarization 
variable simultaneously with the differential cross sec- 
tion involved an interesting tension between the relative 
weight of the x 2 of both observables since the polariza- 
tion is normalized by the differential cross section. The 
results show that in many instances this work (the WP 
OMP) fit the polarization better than it did with the dif- 
ferential cross section. In all three polarization figures 
the major difficulty was again wit h the lighter nuclei, 
specifically the carbon isotopes [4{J, Ill7j | as depicted in 
FigslHH 

The global optical potential of this work is the best in 
reproducing this spin observable over a wide range of tar- 
gets and chains of nuclei. Additional polarization exper- 
imental data is always welcome because it measures the 
interference between the central and spin-orbit terms and 
therefore is a good constraint on their relative strengths. 
The present dataset has a satisfactory amount of these 
reactions at projectile energies lower than 100 MeV but 
additional higher energy polarization experiments would 



be appreciated. 

C. Isoscalar Strength 

To analyze optical potentials it has often been in- 
structive to calculate the volume integrals of the vari- 
ous terms to further illicit theoretical comparison. In 
Fig. [14] the central real and imaginary isoscalar volume 
components at 50 MeV and 150 MeV projectile energy 
are reproduced. There is some ambiguity in determining 
the isoscalar component for the OMPs discussed here so 
the technique used was the average between the proton- 
nucleus and the neutron-nucleus central volume integrals. 
Using this definition the central potential of all tradi- 
tional OMPs can be split into the isoscalar and isovector 
parts respectively 

Vv(E)=V (E)+l(E), (28) 

where Z(E), the isovector component, has its sign de- 
pendent on the isospin of the projectile and Vq(E) is the 
isoscalar term. For Fig. [14] a representative sample of 
nuclei were chosen that are on or near the line of stabil- 
ity, the explicit list of nuclei used in the calculations are 
listed in the figure caption. The green, red, and blue lines 
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FIG. 10: (Color online) Proton-nucleus elastic Rutherf ord reduc ed di fferential cross section experimental data for the target 
nuclei of 12 ' 13 C[1, |U, le O|l 11, 18 ' 20 < 22 O[E M Eil, 2S Si\lO^, 40 < 42 Ca[7l 11, 40 ' 42 < 44 Ar[Il, l|, at varying proton 
laboratory energies. The target labels at the top read left to right correspond to the calculations and experimental data read 
top to bottom for the two leftmost panels. Refer to Fig. [5] for details of the legend for the theoretical calculations. 



are the WP, KD, and MD OMPs volume integrals of the 
isoscalar term (Vo(E) respectively and they have broad 
agreement which is encouraging. All three potentials, fit 
independently to experimental data, developed roughly 
the same strength for the isoscalar central component. 
The only large discrepancy is the MD OMP at low en- 
ergy has a much smaller imaginary component because 
it is missing a surface term. The other two OMPs at 
this energy (KD and WP) have a significant imaginary 
surface term which adds to the strength of the attractive 
central potential. Incidentally if the WP OMP of this 
work is examined closely in Fig. [14] protrusions are seen, 
especially in the imaginary panels. These are due to the 
new magic number term that the potential has included 
(Eqs. I7ll0p . The ramifications of this inclusion are not 
readily apparent in comparisons with the reaction data, 
future studies on the effect of closed shells on optical po- 
tential behavior would be fruitful. 

To elicit further insight the volume integrals created 
from microscopic Watson t — p ( nucleon-nucleon scat- 
tering ma trix folded with target density) potentials [l!8l 
Ill9l . 1 120j ] are plotted in light gray in Fig. [14] The spe- 
cific apparatus was developed, in part, by earlier work 
of on e of the presen t authors (Weppner) and detailed in 
Refs. [HUGH Il23l | In brief, the light gray shaded regions 



use a variety of nucleon-nucleon p oten t ials [124 , Il25l . Il26j | 
folded with theoretical densities [l27l H28j using a vari- 
ety of techniques. These microscopic optical potentials 
have been satisfactory in reproducing experiment at the 
lower energies and better at the higher energies. Rarely 
do the microscopic potentials, which have different fun- 
damental motivations, match the quality of these three 
phcncmological potentials in reproducing these observ- 
ables. Significant systematic differences are recognized 
in the isoscalar volume integrals especially in the light 
nuclei and the imaginary term but there are also regions 
of general agreement. 



In summary, these three OMPs describe the known 
elastic scattering data well and have similar isoscalar 
magnitudes. Detailing their advantages; the KD OMP is 
well suited for neutron projectile projects, the MD OMP, 
in its simplicity, does well with differential cross sections, 
and the WP OMP has strengths in the exotic nuclei and 
in the spin observables. In the next section larger theo- 
retical differences will be proffered, while examining the 
isovector components, that will eventually lead to quite 
disparate scattering observable results. 



18 




30 60 90 30 60 90 30 60 90 120 150 



(deg) (deg) (deg) 

cm. 0/ cm. 0/ cm. °' 

FIG. 1 1: (C olor onli ne) Proton-nu cleus elastic analy zing power spin observable experimental data for the target nuclei of 
28 SiIzl, Eog], 40 Ca[H,S[ii3,[iii], 54 ' 56 Fe[I3,[MIIl[Ili], at varying proton laboratory energies. Refer to Fig.[5]for details 
of the theoretical calculations. 



V. ISOVECTOR TERM ANALYSIS 

The asymmetry isovector term, which measures target 
neutron and proton imbalance, is a point of illustrative 
comparison and will be the focus of this section. This re- 
search attempted to develop an isospin consistent poten- 
tial solely by fitting to elastic experimental reaction data, 
here the success of this effort will be assessed. First the 
historical ansatz of a simple linear N — Z dependence for 
this potential term will be examined, then a Lane anal- 
ysis check of the isospin character of this OMP will be 
proffered. Finally a detailed comparison of the isovector 
volume integral of the three OMPs alongside microscopic 
optical potential results will occur with connections to 
experiment established. 



gin of this N — Z term (following Refs. [HI GM EMI). 
The potential between two nucleons usually contains 
isospin vector components (usually a spin dependent and 
a spin independent piece, for simplification here they are 
combined). Assuming the valid impulse approximation 
at high ene rgies, the tec hniques of Kerman , McManus, 
Thaler [H2 and Watson [Hi ES EM EH can be used 
to construct the nucleon-nucleus optical potential from a 
sum of nucleon-nucleon potentials: 

A 

Vasym ~ ^ VisoN N {^~proj Ti)-> (^^) 

i 

switching over to raising and lowering operators 



A. Validity of N - Z 

Recently the simple linear N — Z dependence of the 
asymmetric potential has been called to question 0. 
All three phcncmological potentials, following tradition, 
have used this standard so it is informative to use micro- 
scopic techniques to develop the optical potential from 
the nucleon-nucleon potential to gain insight into the ori- 



v asym ~ 

A 1 

^ VjsoNN 77 (l~proj-t- ~t~ Tproj — ^"i-f- ^^~proj z ^~i z )• (30) 



If the difference of the asymmetry piece between the pro- 
ton projectile and neutron projectile on the same target 
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FIG. 12: (Color online) P roto n- nucleus elastic analyzing power spin observable experimental data for the target nuclei of 
12 ' 13 C[H,[5l, 16 0[H, 28 Si[lOl, 40 ' 42 < 44 ' 48 Ca[7l[zl, 40 Ar^, 48 < 50 Ti[Il, 50 < 52 ' 54 Cr[zl, at varying proton laboratory energies. 
The target labels at the top read left to right correspond to the calculations and experimental data read top to bottom. Refer 
to Fig. [S] for details of the theoretical calculations. 



is calculated then microscopically it has this format: 

Vasym(,tproj Vasym(tproj ^) 

~v. ( N I \ 

~ VisoNN \ 2 Tproj- ^ ^proj^ ) 
- V isoNN ( CVoja (+ \ ) + T projA-7j))) ■ ( 31 ) 

The first two terms are the simplest inelastic charge- 
exchange term s, often referred to as quasi-elastic charge- 
exchange |l30l ] because they do not involve the direct ex- 
change of nucleons and they are between isobaric analog 
states. The first inelastic term is exclusively non-zero for 
the proton projectile, the second inelastic term is non- 
zero exclusively for the neutron projectile and the last 
two terms are the elastic contributions (the same contri- 
bution occurs for both the proton or neutron projectile). 
A direct N — Z factor is derived directly from the elastic 
scattering component, the inelastic term does not con- 
tain the equivalent proportionality (the proton potential 
contains a non-zero N and the neutron potential contains 
a non-zero negative Z however these potentials are not 
implemented simultaneously). This is the first affirma- 
tion that a linear N — Z term in the optical potential, 



which by definition is both refractive (elastic) and ab- 
sorptive (inelastic), is to first order valid but somewhat 
simplistic. The imbalance of neutrons to protons does 
control the physical mechanisms of the elastic isovector 
piece whereas the number of absolute neutrons (for pro- 
ton scattering) or the number of absolute protons (for 
neutron scattering) are the source of strength for the in- 
elastic component. Although this type of impulse micro- 
scopic approximation does not intrinsically contain mul- 
tiple scattering, correlations, exchange, or coupled chan- 
nels, it is a good approximation at higher energies (it has 
been show n to work adequately at projectile energies of 
150 MeV [120]) and it signifies that beyond first order 
the asymmetry term needs a theoretical re-evaluation. 

Writing the asymmetric potential in macroscopic 
nucleon-nucleus optical potential form is also illuminat- 
ing: 

Vasym — VisoN A\J~proj ' ^~targ) 
— VisoN A\f~projj t -7~t,Q,rg_ ~h Tproj—^~targ-\. H~ ^~proj z ^~targ z ) i 

(32) 

and likewise the difference equation subtracting neutron 
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FIG. 13: (Color online) Prot on- nucleu s elastic anal yzin g power spin observable experimental data for the target nuclei of 
12 ' 13 C[13,[2(l[2l|li, 58 Ni[H,[IM[III|], 59 CoH,m,[l]J|, at var y in g proton laboratory energies. Refer to Fig.[5]for details of 
the theoretical calculations. 



scattering from proton scattering is 

Vasym{tproj Vasym{tproj = ^) 

~ VisoNA^Tproj-Tfarg-i- Tproj-\-Ttarg- TtargJ)- 

(33) 

This form again demonstrates a direct connection to both 
elastic scattering and inelastic charge-exchange. The 
clastic z component piece is proportional to T tar g z — 
N — Z , which confirms the earlier microscopic results. 
Similarly the inelastic piece is again ambiguous. In many 
nuclei the ground state has an isospin vector designa- 
tion such that T tar g z is the maximum value and thus 
\i~targ- \ = l T targ + 1 °c V 'N — Z but wi th u nstable de- 
formed nuclei this is not always the case |l33l | and there- 
fore a general rule about the strength of the inelastic 
piece should be treated w ith apprehension, especially 
at higher energy [l3Ct Il34l |. Both the microscopic and 
macroscopic formulations of the optical potential lead to 
a differentiation between the elastic and inelastic charge- 
exchange compon ent's constant of proportiona lity as seen 
in previous work EH EM EM EM, EM EH . 

All the global optical potentials, perhaps because of 
the linear N — Z term, have trouble giving excellent re- 
sults for the calcium and chromium isotopes (Fig. [5] and 



Fig. UOp . Similarly, as exploration further from the line of 
stability occurs, it is realized that using a stark absolute 
linear term leads to erroneous results. For example using 
a 100 MeV neutron projectile gives physical results when 
scattering from 12 .13,U,15 C using the wp q MP but it 

begins to give negative total cross sections for the same 
scattering observables if the target is 16 C. Obviously all 
elements will return non-physical results when the asym- 
metric term allows for unimpeded linear growth as the 
neutron-proton imbalance gets larger (especially if this 
term is large as is the case for carbon for the WP OMP). 
In general the linear N — Z structure is a good first ap- 
proximation but a better formulation should be devel- 
oped for exotic nuclei optical potentials far from the line 
of stability. 

B. Lane Analysis 

A Lane consistent potential [j fl29l EM fl35l 
guarantees a near equivalent isoscalar and isovector nu- 
clear potential for proton and neutron projectiles at the 
same initial projectile bombarding energies and thus the 
measure of Lane consistency is a good check on the in- 
tegrity of the isovector term. 

As discussed in Sec. IIV C[ the volume potential of 
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FIG. 14: (Color online) A figure of theoretical comparisons. 
Plotted are the isoscalar volume real and imaginary compar- 
isons between the three potentials that are under examination 
at 50 MeV and 150 MeV projectile energy. The light green, 
medium blue, and dark red lines represent the WP, KD, and 
MD OMPs respectively. The light gray shaded area represents 
typical volumes for microscopic t — p potentials (as detailed in 
the text). Each atomic number has a sample nuclei attache 
to it. These nuclei are: 13 C, 14 C, 15 N, 16 N, 17 Q, ls O, 19 F, 



13 C, 14 C, 15 N, 16 N, 17 0, 18 
20 F, 21 Ne, 22 Ne, 23 Na, 24 Na, 25 Mg, 26 Mg, 27 Al, 28 Al, 29 
30 Si, 31 P, 32 P, 33 S, 34 S, 35 C1, 36 C1, 37 C1, 38 Ar, 39 K, 40 Ar, 41 K 

42 n„ 43 



S S CI CI CI ' 

Ca, 43 Ca, 44 Ca,' 45 Sc', 46 Ti, 47 Ti, 48 Ti, 49 Ti, 50 V, 51 V, 52 Cr, 
53 Mn, 54 Fe, 55 Mn, 56 Fe, 57 Co, 58 Ni, 59 Co, 60 Ni. These same 
nuclei will be used in many of the remaining figures. 



all traditional OMPs can be split into the isoscalar and 
isovector parts respectively 



V V (E) = V {E)+1(E), 



(34) 



where 1(E), the isovector component, has its sign depen- 
dent on the isospin of the projectile. Because the effec- 
tive short range projectile kinetic energy of the proton- 
nucleus potential is different than in the corresponding 
neutron-nucleus potential, a Lane consistent potential 
will add a coulomb correction term,A c to the traditional 
proton-nucleus OMP, 



V V {E) =V (E)+1(E)+A C 



(35) 



which adjusts the energy dependent proton-nucleus po- 
tential to account for the lessening of the initial bombard- 
ing energy of the charged projectile as it head s tow ards 
the target due to the long range coulomb field 
be all inclusive this correction should be included for the 
complex volume and spin orbit pieces, traditionally how- 
ever it has only been applie d to the real central term 
(recent exceptions are Ref. [9l. ll34j ). 

Using the notation of Eqs. I5ll27l of this WP potential 
the real central coulomb correction, A c , can easily be 



derived as 

Vv(E)+V s (E) = Vv(E-f cou i)+Vs(E,-f cou i)+A c (E), 

(36) 

where E is the original projectile bombarding energy and 
the volume, surface and coulomb potentials are defined 
using Eqs. OHJ l5l l7|lll This equation leads directly to a 
solution as 

A c = V V (E) + V S (E) - V V (E - f coul ) - V S (E - f coul ), 

(37) 

which is effectively the difference between the original 
and an energy adjusted central volume term. All three of 
the OMPs examined contain a coulomb correction term. 




FIG. 15: (Color online) These are two example depictions at 
65 MeV projectile energy of the functional coulomb correc- 
tion potential for the targets 40 Ca and 56 Fe respectively. The 
solid black line is the magnitude of the exact calculation (the 
coulomb correction is actually negative), for the potential of 
this work, using Eq. 1371 If the coulomb correction term is 
applied to the proton-nucleus optical potential then the po- 
tential can be considered isovector Lane consistent (described 
in the text). The solid green line is the approximate func- 
tion which was determined in this work by fitting the global 
nucleon-nucleus data sets. This approximate form was deter- 
mined by Eq. 1381 The thinner red dashed line is given as an 
example Woods-Saxon functional form for comparison. 

Equation [37] represents the exact definition of the real 
coulomb correction. Historically, starting with Ref. [3J], 
an ansatz was made to use a Woods-Saxon potential form 
with the strength proportional to the average coulomb 
potential in the short range to approximate this coulomb 
correction, both the KD and MD OMP take this ap- 
proach. The WP OMP uses a different approximate 
form, instead of attaching the coulomb correction term 
to the volume term (with Woods-Saxon shape) it is com- 
bined with the coulomb term by having the radius of the 
short term coulomb potential, Eq.[5Sl becomes artificially 
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larger and uncharacteristically energy and atomic num- 
ber dependent. The coulomb correction of this work is 
defined as 



A, 



{Tic) - fcoui{n Co = 1.20 x As/m)) , (38) 



in which a literal addition to the coulomb potential was 
devised by extending IZc beyond its traditional value. 
To extract this coulomb correction to the nuclear term 
we subtracted the known short ranged coulomb potential 
from the full-extended fitted coulomb potential of this 
work as Eq. [581 details. 

This non-traditional approach has some advantages. 
First it disentangles the unknown coulomb correction 
from the likewise unknown nuclear volume term and com- 
bines it to the more apparen t sho rt range coulomb force. 
Second, as developed in Ref. 137| . a Woods-Saxon shape 
is not an excellent functional representation of t he ex act 
coulomb correction result (see Fig. 11 of Ref. [l37l | for 
example). The approximate coulomb correction of this 
work (Eq. [38]) more satisfactorily represents the shape 
of the exact result. Giving an example, in Fig. [TS] are 
depicted two examples at the commonly tested experi- 
mental projectile energy of 65 MeV. The black line rep- 
resents the calculation of the exact result, the green line 
is the function which this research developed in the pro- 
cess of fitting the coulomb radius to the global data set 
using Eq. [35] The thin dashed red line is the Woods- 
Saxon shape normalized to this approximation added as 
a reference. Although the green fit of this research does 
not always match the exact result, its functional shape 
is closer to the exact result specially in the important 
interior region. Another reason to use this approxima- 
tion is that this correction has an operational advantage 
that it can be easy be added to modern optical codes like 
ECIS 15] by changing the value of the coulomb radius. 
The exact result of Eq. [37] which does not have a closed 
analytical form, cannot be added neatly. 

To examine if this inclusion of a coulomb correction 
term brings the WP OMP closer to Lane consistency 
the contour graphs of Figs. [TrjlfTTl are introduced. They 
map over the complete applicable projectile energy and 
atomic number of target phase space for the WP OMP 
and the measure of Lane inconsistency is normalized to 
the moduli of the central potential. The figures have 
ranges between 4% and 10% Lane inconsistency. It is 
noted that of reactions in this phase space the real central 
coulomb correction is at worst about a 12% correction. 
FigHU top panel, shows the level of inconsistency with- 
out a coulomb correction (with an average inconsistency 
of 4.3%). The systematic need for a coulomb correction 
for the present potential is apparent especially at lower 
energies and larger targets, these sensitivities are con- 
gruent with the coulomb correction term developed in 
the KD OMP of Ref. [l[ which increases linearly with Z 
and to good approximation decreases linearly with E. 

The real central coulomb correction is inherently in- 
cluded in this potential, as discussed above, and with it 
the WP potential is substantially more Lane consistent 
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FIG. 16: This is a depiction in projectile energy-target atomic 
number phase space of the real fraction Lane inconsistency 
relative to the magnitude of the total moduli of the central 
potential for sample nuclei. The top panel shows the optical 
potential of this work without a coulomb correction term as- 
suming only 1Z C = 1.20 x As fm. The bottom panel depicts 
the same but with a coulomb correction term assuming a dy- 
namic E and A dependent coulomb radius, 1Z, described in 
the text. The unweighted average for the whole phase space 
of the top panel is 0.043 Lane inconsistency. The unweighted 
average for the whole phase space of the bottom panel is 0.034 
Lane inconsistency. The fraction of 0.04 is assumed here to be 
a pragmatic minimum and a measure of sufficient Lane isovec- 
tor consistency. The sample nuclei for each atomic number 
were near or on the line of stability and are the same as those 
in Fig. EH 



as depicted in Fig. [T|)] bottom panel. Much more of the 
phase space is now below the 4% level and only in a few 
small regions does the coulomb correction addition make 
the potential less Lane consistent (high A and high E 
and low A at midrange energies). The areas where the 
coulomb correction was substandard is where the exper- 
imental data is sparse, especially lacking neutron reac- 
tions. For example, two experimentally popular projec- 
tile bombarding energies are at 65 MeV and 135 MeV. 
Here, with the coulomb correction included (as shown 
in the bottom panel of Fig. [Trj]) . the Lane consistency is 
very good. The energies between these guideposts are 
where the fit has difficulties. The average value of Lane 
consistency with the coulomb correction is 3.4%. In our 
studies developing this analysis it was found that using a 
constant energy independent TZ C between (1.1 x As) fm 
and (1.35 X As) fm to represent the fixed traditional short 
range coulomb radius in Eq. [38] produced the best results 
which is reassuring because it matches the standard short 
range coulomb potential found in the literature. 

The imaginary coulomb correction term has been ig- 
nored in the KD, MD, and this potential. It ha s been ar- 
gued that it can be neglected because it is small |l36l . ll38| 
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FIG. 17: This is a depiction in projectile energy-target atomic 
number phase space of the imaginary fraction Lane inconsis- 
tency relative to the magnitude of the total complex central 
force for sample nuclei which are listed in Fig. 1141 Overall 
the isospin inconsistencies are smaller than its real counter- 
part until E > 140 MeV. The unweighted average for the 
whole phase space is 0.025 Lane inconsistency. The fraction 
of 0.04 is assumed here to be a pragmatic minimum and a 
measure of sufficient Lane isovector consistency. 



The requirement that the optical potential of this work 
be completely Lane consistent is not congruent with the 
goals of this research or the general character of the nu- 
clear force. Cha rge-depen dent (isospin symmetry break- 
ing) discoveries |l39l . Il4fj| ]. the uncertainty in the short 
range coulomb potential, and the approximate functional 
forms applied in the phenemological optical potentials all 
diminish the importance of attaining perfect isospin nu- 
clear symmetry and a Lane potential. The goal of un- 
der 4% inconsistency seems to be a reasonable pragmatic 
minimum giving these conditions. This research aimed 
to be as pure a phenemological fit as possible letting the 
reproduction of experimental data confirm what micro- 
scopic theo ry has shown e lsewhere regarding the coulomb 
correction [1341 1 1371 Il4ll | . The conclusion is that there is 
good agreement when there is a copious amount of exper- 
imental data giving credence to the technique employed. 
C. A comparison 



or because of lack of sufficient experimental data it will be 
ambiguous [l[ . In performing this analysis we found it to 
be small at low energy but unfortunately at high energy 
(greater than 120 MeV projectile energy) it begins to be 
significant in contrast with the real coulomb correction 
term which diminishes as depicted in Fig. [L7j The aver- 
age of the Lane inconsistency over the entire phase space 
normalized to the magnitude of the complex WP central 
potential measures 3.4% in the real term and 2.5% for 
the imaginary term. However the imaginary Lane incon- 
sistency can run higher than 9% for A > 60 and E > 160 
MeV. 

Examining the quality of the results in reproducing ex- 
perimental reactions in Sec. lIVI thcrc is significant error in 
the neutron total cross section calculation at high atomic 
number and high energy as represented in Fig. 3J Since 
this potential fits the proton and neutron observables si- 
multaneously this could be due to the lack of that imag- 
inary coulomb correction term. In future high energy 
optical potential work this imaginary correction should 
not be neglected. 



A comparison of the three global optical potentials 
isovector asymmetry terms, which use the traditional 
N — Z terms, now follows. Here dramatic differences 
between the three formulations can be isolated. 

Starting from the simplest potential, the MD OMP has 
a neutron excess asymmetric term for the real volume and 
the real spin-orbit amplitude only: Q 

N - Z 

Imd = ± — ^ — 16.5 f ws (r,Tli,Ai) 
A dr 

which has no explicit energy dependence and has the 
standard linear term of N ~^ z . The difference between 
the proton projectile and the neutron projectile is sim- 
ply a sign change, the internal geometry parameters have 
no isospin dependence. 

The KD OMP has an explicit asymmetry term only 
for the real volume component and the imaginary surface 
component [l[ 



N — Z 

Xkd = ±—^— (21.0(1 - v 2 (E - E f ) + v 3 (E- E s f - 7.0 x \$-\E - E f ) 3 )f ws (r,TZ lv ,A lv )) 

Ti±A d ^^-{^ f J E ~ E 2 f ] J 2 ex V (-d 2 (E - E f ))^-(f ws (r,n d ,A d )), (40) 
A (E — Ef) + «3 



where v%, V3 , di and ds are functions which depend on 
the nucleon number of the target, and projectile energy 
and the isospin character of the projectile. Likewise Ef 
represents the Ferm i en ergy of the target, extracted from 
mass excess values jl42| , and is dependent on the isospin 



of the projectile. The internal geometry parameters have 
no explicit isospin dependence. Because of the separate 
functional dependence on the projectile the asymmetry 
term, Ir- jj> of Eq. 1401 is not exactly linear and the isospin 
flip in the projectile is not simply a sign change as with 
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the Imd term given by Eq. [39l 

The optical potential of this work (WP OMP), out- 
lined in Sec. Ill Al has vector isospin asymmetry (N — 
Z) in five major terms: real and imaginary volume, 
imaginary surface, and real and imaginary spin orbit 
( , Eas. l6l9ll3lT5TT7ll : 

1 WP = ±(N - Z){V Vi + iW Vi )fws(r,Ki,Ai) 
Ti4(N - Z)A s W St ^fwsir, n s ,A s ) 

T(N - Z){V SOi +iW S o i )^-fws(r,n S o,Aso) (41) 

dr 

where Vy 4 , Wyj , , Vso t , and WsOi are separable poly- 
nomial functions in terms of projectile energy and nu- 
cleon number (Eqs. [5"1 [2~T)) . There is no projectile isospin 
dependence within the polynomials. This potential also 



Jv v /A = 
Jw v /A = 
Jws/A = 



These equations will be used to calculate the differences 
in volume integrals between a proton projectile and a 
neutron projectile acting upon the same target nucleus 
(the isospin asymmetry). The difference is emphasized, 
as it was in the linear study of the N — Z term, because 
this is what characterizes the isovector asymmetry term 
from the rest of the dominating isoscalar nuclear poten- 
tial; it does not disappear upon subtraction of this isospin 
flip of the projectile. In the KD potential, where there 
are two separate potentials with different functions to dif- 
ferentiate proton and neutron scattering, this definition 
of asymmetry is somewhat ambiguous, it is much more 
enhanced than simply the N — Z term. So in this com- 
parison we will define isovector to be half the difference 
between the proton-nucleus and neutron-nucleus poten- 
tial. This will include the addition of the nuclear coulomb 
correction term in the proton potential for all three of the 
OMPs. This seems to be the best workable definition for 
isovector for these three potentials especially in tandem 
comparisons with the microscopic potential that follow. 

In Fig. [18] the two panels contain results of the real 
volume, surface, and coulomb correction integral differ- 
ence at 50 MeV (left panel) and 150 MeV (right panel) 
energy between a proton and a neutron projectile for the 
three global optical potentials (approximately twice the 
isovector strength). This volume integral difference is a 
direct measure of the change in strength that the three 
potentials have for the dynamics of a isospin flip of the 



enforces the separability of E and A by using (N—Z) and 
not N ^ Z ■ There are some explicit asymmetry terms in 
the Av (Eq.[20|) and Hso (Eq. 123)1 geometry terms which 
also lead to non-symmetric neutron excess terms. This 
optical potential has therefore attempted to fit the imag- 
inary volume, surface, and spin orbit asymmetry terms, 
the other contemporary global OMPs discussed within 
have set their explicit imaginary asymmetric terms to 
zero. 

To explore the differences between the three OMPs 
further let us again examine volume integrals which have 
been illustrative in the past to help describe Gamow- 
Teller and Fermi charge-exchange transitions |l43| and 
to effectively compare the strength of disparate shaped 
optical potentials. Explicitly calculated are, using the 
notation of Eq. [2) the following integrals: 



(42) 
(43) 
(44) 



projectile on the same N — Z ^ target. For every 
nucleon number (13 < A < 60) there was chosen a rep- 
resentative target that is stable or close to the line of 
stability which are listed for Fig. [T?J The calculation is 
explicitly of 

(Re J/A) iao = (Jv v /A) p 

roton 

(JvJA) 

neutron 

(45) 

using the Eq. [15] and the integral of the coulomb cor- 
rection term, Joeita a , which for the WP OMP has the 
analytical form 

2 / 2 \ 

Ja c = --irhceZ [pl c ) 2 - (1.20 x Aifm) j , (46) 

and for the KD and MD potential the coulomb correc- 
tion integral calculation is similar to Eq. 05] because of 
its Woods-Saxon functional form. Figure [TSl shows some 
agreement at the low energy (for 28 < A < 60) but it is 
more disparate at the high energy. Likewise at A < 20 
there is sharp disagreement between the two applicable 
potentials (MD and WP) for these light nuclei. The fine 
features are telling also, the wildly oscillatory behavior 
of the MD and WP integrals are indicative of the size of 
the neutron excess (this is to be expected because they 
are similar to Gamow- Teller and Fermi sum rules which 
are proportional to N — Z [1431 ]), in the KD integrals 
this behavior has been quenched because the neutron and 
proton potentials lack the requisite similarity. 



~J r 2 (y v (E,A,N,Z,V,MN))f W s(r,K v ,Av) dr 
r 2 (W v [E, A, N,Z,V, MN)j f ws (r, K V ,A V ) dr 
J r 2 (w D (E,A,N,Z,V))-^fws(r,n Sl As) dr. 
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FIG. 18: (Color online) This is a calculation of the difference 
in the real component of the nuclear central potential inte- 
grals (volume and coulomb correction) between a proton and 
neutron projectile (about twice the actual isovector volume 
integral) for a representative selection of nuclei using three 
different global optical potentials. The nuclei were chosen to 
have nonzero N — Z terms and to be close or on the line of 
stability. The nuclei are the same as used in Fig. 1141 The dark 
blue semi-filled circular points are the calculation using the 
KD optical potential [JJ], the red unfilled circular points use 
the MD optical potential Q] and the light filled green circu- 
lar data points are calculated using the WP potential of this 
work. The lines connecting the data points are there to make 
it easier for the eye to follow. The left panel used a 50 MeV 
nucleon projectile and the right panel uses a 150 MeV pro- 
jectile where the differences in the three calculations become 
more pronounced. The light gray shaded area is the range of 
values for microscopic potentials which are described in the 
text in Sec. UVCl 



for the source of the magnitude shifts are not substan- 
tially isovector in origin but are almost entirely due to 
the coulomb correction which is very much reduced in 
the WP and KD optical potentials at high energies but 
is energy independent in the MD optical potential. 
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FIG. 19: (Color online) This is a calculation of the differ- 
ence in the imaginary component of the central volume inte- 
grals between a proton and neutron projectile (about twice 
the isovector strength) for a representative selection of nuclei 
using the WP and KD global optical potential. The nuclei 
were chosen to have nonzero N — Z terms and to be close or 
on the line of stability, the specific nuclei are listed in Fig. Q3] 
The light green fill represents the range possible for nucleon 
projectile energies between 50 MeV and 150 MeV for the WP 
global optical potential of this work. The KD global opti- 
cal potential (dark blue fill) is smaller and has a diminished 
range but its values are comparable to its real counterpart. 
The light gray shaded area is the range of values for micro- 
scopic potentials which are described in the text in Sec. IIV Cl 



Figure [TH] also has a light gray shaded area depicting 
a range of microscopic optical potentials isovector vol- 
ume integrals. These microscopic potentials have been 
described in Sec. IIV CI and their isoscalar volume inte- 
grals depicted in Fig. [14] Overall the microscopic optical 
potentials have a systematically lower strength than the 
three OMPs, although it is encouraging to see that at 
the lower 50 MeV projectile energy there is a general 
agreement between all the potentials for target nucleons 
with A > 28. Some features of the microscopic poten- 
tial isovector character are that it has a much smaller 
width at higher energies because the impulse approxima- 
tion, where a free density-independent t-matrix is closer 
to being realized. Likewise the microscopic potential also 
shows an oscillatory N — Z behavior in the isovector vol- 
ume strength, for example the isovector character of 14 C 
is almost twice as big as 13 C because it has twice the neu- 
tron excess. At high energies the KD and WP isovector 
elements shrink and mimic the microscopic result bet- 
ter than the MD OMP which grows; this is deceiving, 



The calculation of the explicit imaginary volume asym- 
metric potential is a distinct attribute of the WP optical 
potential described by this work. Calculating the isospin 
difference of the imaginary term for this potential is ex- 
ecuted by: 

(ImJ/A)i s0 — (Jw v /A + Jws/A) proton 

- {Jw v /A+Jw s /A) (47) 

where both the volume and surface terms are included 
as described by Eqs. |4~2"R41 In contrast the KD OMP 
only has an implicit small imaginary asymmetric sur- 
face term but in subtracting a neutron projectile reac- 
tion from a proton projectile reaction off the same target 
nucleus there is a significant imaginary difference coming 
from the disparity in isospin dependent volume strengths 
and Fermi energies. The MD potential has an imaginary 
isovector strength of zero. In Fig. [TH] a plot of these 
imaginary volume integrals are depicted at a range be- 
tween 50 MeV and 150 MeV projectile energies. The re- 



20 



suits are striking, this imaginary volume integral is quite 
large, often the same size or larger than the companion 
real piece. Having an imaginary asymmetry term, as the 
WP and KD optic potentials have, allows for a more nu- 
anced optical potential which contains mechanisms for 
charge-exchange resonances |l43l | , important at these en- 
ergies. Again, for further comparison, the microscopic 
optical potential range for the imaginary isovector ele- 
ment is depicted in light gray. It is systematically lower 
than the phenemological optical potentials but is also is 
often much stronger than its real counterpart, especially 
at higher energies. The microscopic potential integral 
still carries large remnants of the N — Z character in the 
results as does the WP OMP. The KD OMP is greatly 
reduced reminiscent of the real isovector volume element 
of Fig. 1181 Incidentally, if a realistic imaginary coulomb 
correction is added it would give a small positive volume 
element differential to the results depicted in Fig.lT9lthus 
leading to better agreement between the phenemological 
and microscopic potentials for this term. 

A comparison of these asymmetric isovector volume 
integrals differences can me made to experiment. It is 
well understood that the asymmetry volume integral is 
tied to Gamow- Teller transitions in charge- exchange re- 
actions which dominate at high energies 143]. The vol- 
ume integral equations of Eqs. l42ll47l are equivalent to do- 
ing a Fourier transform to momentum transfer (q) space 
set to zero. In this forward scattering case the inelastic 
aspect of the asymmetry terms dominates. This is be- 
cause in the coulomb distorted wave basis the long range 
coulomb potential, which dictates extremely forward an- 
gle proton scattering, is external to the optical potential 
so the magnitude of the elastic scattering component is 
zero, likewise the neutron scattering elastic forward am- 
plitude is zero. This extreme scattering has been recog- 
nized as b eing imp ortant in determining Gamow- Teller 
strengths [l44L 1 1451 ] . If the difference is taken between 
the energy needed to initiate a proton-neutron charge- 
exchange and a neutron-proton charge-exchange for the 
same target nucleus (the difference in Q values derived 
from the mass excess values for the resultant final nuclei) 
this can be compared to the difference in the volume in- 
tegrals which give the minimum energy, zero momentum 
transfer, nuclear density difference. 

The volume integral for the coulomb potential is infi- 
nite but it can easily be defined within the traditional 
nuclear radius, using the short range expression, used 
by all three OMPs under examination, which states that 
when r < He- 

71 A~ 

Jc/A= ^lJ V fc,ui(r,K c ,A,N,Z)r 2 dr. (48) 

In the short range the coulomb potential interferes with 
the nuclear force and contributes to the overall strength 
of the asymmetric volume element. The KD and MD 
optical potentials also include a separate coulomb cor- 
rection term while the potential of this work (WP) does 
not contain a separate coulomb correction term but is 



part of the original short-range coulomb term as detailed 
in Sec. El These nuclear corrections are significant in the 
short range and must be included when illustrating short 
range differences between neutron and proton scattering. 

To calculate the difference per nucleon in potential 
energy density at zero momentum transfer in the short 
range five terms are calculated 

Jiso \R& J /A -(- Jcoul.corr. /A^ so 

+ (Im J/A) iso 

+ J c /A-AJ sr , (49) 

where the first line is the real nuclear potential energy 
density difference with the coulomb correction included 
(Fig. [T5]). The second line is the imaginary nuclear po- 
tential energy density difference (Fig. [THl The last line 
is the short range coulomb potential volume (Eq. [48]) , up 
to the nuclear radius 7?.y, and a zero point correction, 
AJ sr . This correction is needed because of the weak- 
ened short ranged coulomb volume that is used by all 
three potentials as dictated by Eq. [4] Since the coulomb 
distorted wave basis sets zero nuclear scattering as de- 
fined by the traditional long range coulomb form (Eq. [3]) 
and the short range volume element form actually used 
is 20% weaker then the traditional coulomb potential, 
the true zero-scattering point has been shifted in the dis- 
torted wave coulomb basis of neutron and proton scat- 
tering within the short range and this modification has 
to be normalized accordingly. Thus the correction is 

AJ sr - (0.2)4 / — r 2 dr, (50) 

A Jo r 

which is applied to all the potentials consistently. 

In Fig. [201 and Fig. [5T] there are depicted Q value 
(mass excess) differences for charge-exchange reactions 
(Q(n,p) — Q{p, n)) in black along with normalized volume 
integrals for the MD, W P, K D, and a typical microscopic 
t — p optical potentials jl22| for the forty-eight different 
sample target nuclei at a projectile energy of 150 MeV. 
The normalized optical potentials were calculated by us- 
ing Eq. 1491 and then multiplying that result by a constant 
density. The functional form results of the optical poten- 
tial volume integrals are in remarkable agreement with 
the experimental Q value differences. Three of the four 
optical potentials mimic the shape and structure of the 
experimental mass excesses, following the global mini- 
mum at 13 C, local maxima and minima at 39 K and 40 K 
(dictated by factors of N — Z), and the general trends 
as nucleon number increases. A natural consequence of 
fitting elastic scattering data with a proton and neutron 
inclusive optical potential is the ability to mimic the Q 
values of the charge-exchange reactions at high energies 
given a simplistic constant density normalization factor 
even without this data being used as a fitting constraint. 
This structure is remarkably also built into the the micro- 
scopic optical potential which is a sum of the two-body 
nucleon- nucleon potentials jl22| . The only optical poten- 
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FIG. 20: (Color online) This plot depicts the total normal- 
ized energy difference (Ji ao pnuc) for the proton and neutron 
projectiles on a representative selection of nuclei using two dif- 
ferent global optical potentials using potentials at projectile 
energy of 150 MeV. The nuclei were chosen to have nonzero 
TV — Z terms and to be close or on the line of stability, the 
specific nuclei are listed in Fig. 1141 These calculations are 
compared to the difference in experimental mass excess ener- 
gies for the charge- exchange reactions, Q(n,p) — Q(p,n) from 
the same targe ts, t he experimental data (dark black circles) 
are from Ref. [146| . The MD calculations, the unfilled red 
circles, have multiplied the results of Eq. [49] by .38 /m -3 , 
the WP calculations, the light filled green circles, has been 
multiplied by .16 /m -3 . The experimental mass excess ener- 
gies shown here were not used to constrain the OMPs during 
the fitting procedure. Two other potential calculations along 
with experiment are shown in Fig. [21] 



tial that does poorly is the KD potential which was de- 
veloped using separate potentials for proton and neutron 
scattering. The general conclusion is that by simultane- 
ously fitting neutron and proton projectiles (or neutron- 
proton and proton-proton phase shifts in the case of the 
microscopic potentials) one automatically develops a po- 
tential which has the correct target dependant relative 
moduli strengths for the forward scattering volume inte- 
grals at high energy which are tied directly to both the 
mass difference and Gamow- Teller matrix elements. 

The general systemic target shape is correct in 
Figs. 1201211 for the MD, WP, and t - p potentials but 
the magnitude differences are substantial. The constant 
density factor is 1.14 fm -3 for the microscopic potential, 
.38 fm -3 for the MD calculation (the isovector volume 
strength is three times the size of the microscopic calcu- 
lation) and the WP OMP has .16 fm~ 3 for the constant 
density (the isovector volume strength is seven times the 
size of the microscopic) . Additionally the MD optical po- 
tential has all its strength in the real component where 
the WP and microscopic potentials have isovector charac- 
ter distributed in both the real and imaginary term. Oth- 
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FIG. 21: (Color online) This plot depicts the total normal- 
ized energy difference (Ji SO pnuc) for the proton and neutron 
projectiles on a representative selection of nuclei using two dif- 
ferent global optical potentials using potentials at projectile 
energy of 150 MeV. The nuclei were chosen to have nonzero 
N — Z terms and to be close or on the line of stability, the spe- 
cific nuclei are listed in Fig. 1141 These calculations are com- 
pared to the difference in experimental mass excess energies 
for the charge-exchange reactions, Q(n,p) — Q(p,n) from the 
same targ ets, t he experimental data (dark black circles) are 
from Ref. [l46j . The KD calculations, the filled blue circles, 
have multiplied the results of Eq.[49]by .17 /m -3 , the light er 
gray circles (following the technique described in Ref. |122| ]). 
has been multiplied by 1.14 /m -3 . The experimental mass 
excess energies shown here were not used to constrain the 
OMPs during the fitting procedure. 

ers have used charge-exchange differential cross-sections 
to constrain the nucleon-nucleus optical potential asym- 
metric isovector term. Future work on this OMP could 
take t hat direct ion following the procedures as outlined in 
Refs. [l34l Il47l | which would help constrain these isovec- 
tor magnitudes further. 

The ramifications of the large asymmetric potential 
differences can be ascertained also with reactions that 
the global optical potentials are currently fitted to. Ta- 
ble IIVI examines the differences between the 40 Ar(n,*) 
and 40 Ca(n,*) cross sections with a 50 MeV and 150 
MeV neutron projectile (calcium subtracted from argon). 
The large differences generated in the elastic and inelastic 
cross section predictions are highly dependent on the con- 
trasts in the antisymmetric isovector term between the 
three potentials. This reaction was chosen because both 
A — 40 nuclei are stable and there are dramatic differ- 
ences between the potentials at A — 40 in the asymmetric 
term. Using neutrons as projectiles also allows an igno- 
rance of coulomb effects and coulomb correction terms. 
At the 50 MeV projectile energy the isovector differences 
are large however the calculations fall within on e or two 
standards deviation of experimental error [9^, Il48| | At 
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Energy 


Model 


AReJ 


AlmJ 


A<TeZiist 


A(7react 


Acr tot 


MeV 




MeV fm 3 


MeV fm 3 


mb 


mb 


mb 




KD 


+ 1.8 


+0.4 


+ 14 


-28 


-14 


50 


MD 


+ 14.5 


0.0 


-36 


-9 


-45 


MeV 


WP 


-10.8 


+8.2 


+68 


-51 


+17 




t-p 


+14.0 


+7.5 


-18 


-22 


-40 




KD 


+2.6 


+0.5 


-19 


-4 


-23 


150 


MD 


+18.2 


0.0 


-106 


-5 


-111 


MeV 


WP 


+3.1 


+16.8 


+188 


-287 


-99 




t-p 


+2.1 


+5.8 


-6 


-15 


-21 



TABLE IV: This is a comparison of the strength of the asymmetry term at A = 40. The calculation is the difference between 
neutron-nucleus total cross section calculations from scattering off of 40 Ar and 40 Ca at 50 MeV and 150 MeV projectile energies 
{p Ar — o"Ca). The differences between the three optical potentials at high energies are mostly due to the isovector asymmetric 
term and pronounced and easily separable. At lower energies the differences fall within experimental deviation. 



the higher projectile energy of 150 MeV, where it can be 
assumed that the born and impulse approximations have 
validity, the effects of correlations, coupled channels and 
potential distortions are at a minimum. The extreme dra- 
matic reduction in the inelastic cross section and increase 
in the elastic cross section predicted by the WP OMP dif- 
fers wildly with the predicted large decrease in the elastic 
cross section seen by the MD potential. Likewise the KD 
and microscopic potentials offers a third very realistic 
conclusion; that the differences between 40 Ar(n,*) and 
40 Ca(n,*) cross sections with a 150 MeV neutron pro- 
jectile are relatively minor. There is at present no high 
energy 40 Ar(n,*) data to confer which optical potential is 
closest to experiment. Experiments calculating the 40 Ar 
cross sections would be a welcome addition as well as 
would other high energy neutron rea ction data like re- 
cent experiments detailed in Refs. [Ill Ill0l | which, if use 
N Z target nuclei, are sensitive to the isovector term. 



VI. CONCLUSION 



also ascertained that a benefit of fitting the proton and 
neutron observables simultaneously is the ability to ac- 
curately define the structure of the mass excesses for the 
charge-exchange reactions. 

There also was a critical examination of the validity of 
the traditional linear proportionality of N — Z applied in 
the asymmetry term. An examination of the calculation 
of the calcium and chromium experimental observables 
produced at 65 MeV projectile energy show a continuing 
failure to produce great fits along the isotopic chain, in- 
dicative of a breakdown of this linear assymetry anasatz. 
Further systematic experimental studies of chains of iso- 
topes would help ascertain the direction of further theo- 
retical studies. 

Future work could extend this potential to heavier tar- 
gets, test other forms of asymmetry potentials, examine 
spin-spin terms, use charge-exchange information to con- 
strain data and determine algorithms to better weight 
the scattering observables. To better constrain the terms 
more elastic scattering data from traditional and exotic 
nuclei data is needed especially at the higher energies. 



The motivations for this work were to construct a 
phenemological nucleon-nucleus global optical potential 
that is suited for a wide range of nuclei targets and pro- 
jectile energies which are within capacity of the exotic 
beam accelerators presently running and under develop- 
ment. We have succeeded in creating one isospin depen- 
dent potential which fits target nuclei 12 < A < 60 and a 
projectile energy of 30 MeV< E < 160 MeV. It compares 
well with two other recent global optical potentials [H, Q , 
its advantages are that it is one continuous optical poten- 
tial which is designed to do systematic studies on mirror 
nuclei and chains of isotopes (a observable calculator has 
also been made available to researchers to quickly use 
this potential for their own research [IH). We have also 
included an imaginary asymmetry term which is miss- 
ing from other recent global optical potentials and have 
given comparative analysis on how the asymmetric po- 
tential terms in the three potentials dramatically differ 
from each other, microscopic potentials and the ramifi- 
cations to experimentally testable observables. It was 
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